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ABSTRACT 




We have calculated the two-point correlation functions in redshift space, £(c, tt), for galax- 
ies of different spectral types in the 2dF Galaxy Redshift Survey. Using these correlation 
functions we are able to estimate values of the linear redshift-space distortion parameter, 
(3 = fl^f /b, the pairwise velocity dispersion, a, and the real-space correlation function, £(r), 
for galaxies with both relatively low star-formation rates (for which the present rate of star 
formation is less than 10% of its past averaged value) and galaxies with higher current star- 
formation activity. At small separations, the real-space clustering of passive galaxies is very 
much stronger than that of the more actively star-forming galaxies; the correlation-function 
slopes are respectively 1.93 and 1.50, and the relative bias between the two classes is a de- 
clining function of radius. On scales larger than 10ft. -1 Mpc there is evidence that the relative 
bias tends to a constant, 6 paS sive /^active — 1- This result is consistent with the similar degrees 
of redshift-space distortions seen in the correlation functions of the two classes - the contours 
of £(er, tt) require /3 a ctivc = 0.49 ± 0.13, and /3 paS sivo = 0.48 ± 0.14. The pairwise velocity 
dispersion is highly correlated with (3. However, despite this a significant difference is seen 
between the two classes. Over the range 8 — 20 h^ 1 Mpc, the pairwise velocity dispersion 
has mean values 416 ± 76 kms -1 and 612 ± 92 kms -1 for the active and passive galaxy 
samples respectively. This is consistent with the expectation from morphological segregation, 
in which passively evolving galaxies preferentially inhabit the cores of high-mass virialised 
regions. 



Key words: galaxies: statistics, distances and redshifts - large scale structure of the Universe 
- cosmological parameters - surveys 



2 D.S. Madgwick et al. 



1 INTRODUCTION 



2 THE 2dFGRS DATA 



It is now well established that the clustering of galaxies at low red- 
shift depends on a variety of factors. Two of the most prominent 
of these, which have been discussed extensively in the literature, 
are luminosity (see e.g. Norberg et al. 2001 and references therein) 
and galaxy type (e.g. Davis & Geller 1976; Dressier 1980; Lahav, 
Nemiroff & Piran 1990; Loveday et al. 1995; Hermit et al. 1996; 
Loveday, Tresse & Maddox 1999; Norberg et al. 2002a). It is the 
latter of these which we wish to address in this paper, by making 
use of the 221,000 galaxies observed in the completed 2dF Galaxy 
Redshift Survey (2dFGRS; Colless et al. 2001). 

Previous analyses of the clustering of galaxies as a function 
of morphological type have revealed that early-type galaxies are 
generally more strongly clustered than their late-type counterparts; 
their correlation function amplitudes being up to several times 
greater (e.g. Hermit et al. 1996). These results are also found to 
hold true if one separates galaxies by colour (Zehavi et al. 2001) 
or spectral type (Loveday, Tresse & Maddox 1999), both of which 
are intimately related to the galaxy morphology (see e.g. Kenni- 
cutt 1992). 

The existence of this distinction in the clustering behaviour of 
different types of galaxies is to be expected if one considers galax- 
ies to be biased tracers of the underlying mass distribution, since 
the amount of biasing should be related to a galaxy's mass and for- 
mation history. However, the fact that the most recent analyses of 
the total galaxy population have revealed that galaxies are not on 
average strongly biased tracers of mass on large scales (Lahav et 
al. 2002; Verde et al. 2002) makes this behaviour even more inter- 
esting, and puts some degree of perspective on these results. 

In this paper we attempt to make the most accurate measure- 
ments of this distinctive clustering behaviour, by calculating the 
two-dimensional correlation function, £(er, 7r), for the most quies- 
cent and star-forming galaxies in our sample separately, where a 
is the galaxy separation perpendicular to the line-of-sight and n 
parallel to the line-of-sight. This simple statistic allows us to eas- 
ily visualise and quantify the variation in clustering properties on 
a variety of scales, picking out for example the distinctive 'finger- 
of-God' effect due to peculiar velocity dispersions in virialised re- 
gions, as well as the large scale flattening due to coherent inflows 
of galaxies towards over-dense regions. 

By contrasting these observed effects we can gain significant 
insights into the properties of the galaxy population, particularly 
when these results are set against simple analytic models. For ex- 
ample, we can use the large-scale inflows to constrain the quan- 
tity j3 = Q. 6 /b, and the small-scale 'finger-of-God' distortions 
to constrain the distribution of galaxy peculiar velocities f(v) si- 
multaneously. Such an analysis has already been performed us- 
ing an earlier subset of the 2dFGRS by Peacock et al. (2001), and 
an updated version of that analysis is presented in Hawkins et al. 
(2003). This paper extends their analyses by incorporating the spec- 
tral classification of 2dFGRS galaxies presented in Madgwick et al. 
(2002a). 

The outline of this paper is as follows. In Section 2 we briefly 
describe the 2dFGRS data-set and the spectral classification we are 
adopting. In Section 3 we then outline the methods we use for esti- 
mating the correlation function and the models we use when mak- 
ing fits to this function. The results of our parameter fits are pre- 
sented in Section 4 and are compared to previous results in Sec- 
tion 5. In Section 6 we conclude this paper with a discussion of our 
results. 



The data-set used in this analysis consists of a subset of that pre- 
sented in Hawkins et al. (2003) - including only those galaxies 
with spectral types (Section l2.lt . which lie in the redshift interval 
0.01 < z < 0.15. Again we restrict ourselves to only considering 
the most complete sectors of the survey, for which > 70% of the 
galaxies have successfully received redshifts. This leaves us with 
96 791 galaxies for use in the analysis presented here. Further de- 
tails of the data-set are presented in Hawkins et al. (2003). 



2.1 Spectral types 

We adopt here the spectral classification developed for the 2dF- 
GRS in Madgwick et al. (2002a). This classification, rj, is derived 
from a principal components analysis (PCA) of the galaxy spectra, 
and provides a continuous parameterisation of the spectral type of 
a galaxy based upon the strength of nebular emission present in its 
rest-frame optical spectrum. It is found that rj correlates relatively 
well with galaxy B-band morphology (Madgwick 2002). How- 
ever, the most natural interpretation of rj is in terms of the relative 
amount of star formation occurring in each galaxy, parameterised 
in terms of the Scalo birthrate parameter, bscaio (Scalo 1986), 



bScalo 



SFR 



present 



(SFR) 



(1) 



past 



as demonstrated in Madgwick et al. (2002b). 

Although 77 is a continuous variable we find it convenient to 
divide our sample of galaxies at r\ = —1.4. It is found that this cut 
of rj = —1.4 corresponds to approximately bscaio = 0.1 (i.e. the 
current rate of star formation is 10% of its past averaged value). 

The cut in rj we have adopted is the same as that used to dis- 
tinguish the so-called 'Type 1' galaxies used in our calculation of 
the galaxy bj luminosity functions (Madgwick et al. 2002a). In that 
paper two more cuts were made at rj — 1.1 and rj = 3.5, which 
we have not adopted in this analysis. It is found that using only two 
spectral types instead of four greatly increases the accuracy of our 
analysis, while the clustering properties of the most actively star- 
forming galaxies are found to be very similar (Section 13.21 . The 
clustering with spectral type in the 2dFGRS has previously been 
considered by Norberg et al. (2002a); however the present paper 
extends this analysis by considering the full magnitude limited sur- 
vey. 

In the analysis that follows the two samples constructed by 
dividing at 77 = —1.4 will be referred to as the relatively passive 
and active star-forming galaxy samples. These two samples consist 
of a total of 36318 and 60473 galaxies respectively. 



3 THE TWO-POINT CORRELATION FUNCTION 

The correlation function, £, is measured by comparing the actual 
galaxy distribution with a catalogue of randomly distributed galax- 
ies. These randomly distributed galaxies are subject to the same 
redshift and mask constraints as the real data. £(<r, n) is estimated 
by counting the pairs of galaxies in bins of separation along the 
line-of-sight, 7r, and across the line-of-sight, a, using the following 
estimator, 



fa,*) 



(DP) - 2{DR) - (RR) 
(RR) 



(2) 
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Figure 1. The £(cr, 7r) grids for our different spectral types: (a) passive, (b) active and (c) full samples. The contour levels are at f = 4.0, 2.0, 1.0, 0.5, 0.2 and 
0.1. 



from Landy & Szalay (1993). In this equation (DD) is the 
weighted number of galaxy-galaxy pairs with particular (a,n) sep- 
aration, (RR) the number of random-random pairs and (DR) the 
number of data-random pairs. The normalisation adopted is that 
the sum of weights for the real galaxy catalogue should match that 
of the random catalogue as a function of scale. As in Hawkins 
et al. (2003) we adopt the J3 weighting scheme to minimise the 
variance in the estimated correlation function (Peebles 1980). We 
have also estimated the correlation functions using the estimator of 
Hamilton (1993), however because these give an essentially identi- 
cal estimate for £ (well within the statistical uncertainties) we only 
present the results from the Landy & Szalay estimator in this paper. 

The random catalogue is constructed by generating random 
positions on the sky and modulating the surface density of these 
points by the completeness variations of the 2dFGRS. Note that, in 
contrast to Hawkins et al. (2003), this completeness now also in- 
cludes the that introduced by the spectral classification for galaxies 
with z < 0.15 (Norberg 2001). The redshift distribution is then 
drawn from the selection function of each type as calculated from 
the 2dFGRS luminosity functions - which allow us to naturally 
incorporate the varying magnitude limit of the survey. These lumi- 
nosity functions have been calculated as in Madgwick et al. (2002a) 
for the data used in this analysis in both the NGP and SGP regions 
separately for each spectral type. 

Due to the design of the 2dF instrument, fibres cannot be 
placed closer than approximately 30 arcsec (Lewis et al. 2002). In 
Hawkins et al. (2003) the effects of these so-called 'fibre collisions' 
were taken into account in the estimation of the correlation func- 
tions by comparing the angular correlation functions of the parent 
and redshift catalogues. It was found that the effect was signifi- 
cant for separations < 0.2 ft -1 Mpc. This cannot be done for the 
present analysis because of the spectral type selection and so we 
ignore all separations < 0.2 /i -1 Mpc in the fitting process. 

The resulting estimates of £(<r, it) are presented in Fig.Qand 
clear differences are immediately visible. The rest of this paper is 
spent quantifying these differences. 



3.1 De-correlating error bars 

Many of the subsequent sections of this paper will involve attempt- 
ing to fit a parametric form to £(<r, n) or f (r). Because the individ- 
ual points we estimate for each of these quantities and their associ- 



ated error bars are not independent (a single galaxy can contribute 
to correlations over all scales), a simple % 2 fit between the observed 
and model correlation functions may not yield the most accurate 
parameter fit. For this reason we must carefully account for these 
correlations in each of our subsequent fits. 

When fitting the model correlation function to the observed f , 
we are interested in minimising the residual between the two. For 
this reason we make a simple change of variables to define, 

(«*) - 



A( Si ) = 



(3) 



Where here £(s) is a given realisation of the correlation function 
we are estimating, (£(si)) is the mean value of our ensemble of 
correlation functions at separation Si, and <T((si) is the standard 
deviation of the estimates of the correlation function at this same 
separation, as deduced from the bootstrap analysis described below. 
We then construct the covariance matrix, 



( 1 



(A( Si )A( 8j )) 



(4) 



in terms of these variables. 

The best-fitting model correlation function £ modcl (s) can be 
found through minimising the residual between it and the observed 
correlation function £ obs (s) in terms of the x 2 difference between 
the two. The residual between the models and observations is de- 
fined by, 



v ( res )/ 



/ Aobs / „ \ /^model/„ \\ 
cr e (si) 



(5) 



in which case the \ 2 can be found from, 



x - = (A ylua >) 1 C- 1 A {r ' lsl . (6) 

Where A' rcs ' is the vector of elements given in Eqn.|5| The above 
equations can then easily be generalised for the gridded £(cr, n). 
Because the points in the £ (a, tv) grid are highly correlated there 
are in fact very few independent components in the observed cor- 
relation function. For this reason it is also possible to extend this 
analysis as demonstrated by Porciani & Giavalisco (2001) by in- 
stead first performing a principal component analysis using our es- 
timated covariance matrix, however we have not found this step to 
be necessary for this analysis. 

One would require a set of independent realisations in order 
to determine unbiased estimates of <r^ and the data covariance ma- 
trix C. However, this is of course not possible since we have at 



1 X (res) 
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our disposal only one realisation of the Universe. In most cases 
a good determination of this covariance can nonetheless be deter- 
mined from mock galaxy simulations (e.g. Cole et al. 1998). These 
simulations are very useful in that they represent the total expected 
cosmic variance. 

The complication we encounter as opposed to other analyses 
(e.g. Hawkins et al. 2003) is that the simulations available to us 
cannot adequately account for the variation in galaxy clustering for 
different types of galaxies. For this reason these simulations can 
only provide us with rough estimates as to the magnitude of the 
cosmic variance - which we must somehow scale to correspond to 
our results. 

Another method for determining the expected covariance of 
our data is to use a bootstrap re-sampling of the data-set (Ling, 
Barrow & Frenk 1986). The most important assumption in a boot- 
strap estimation of the covariance matrix is that each of the data- 
points sampled must be independent. This is not true of the galaxy 
distribution itself, but if we divide our survey area into a selection 
of contiguous regions and re-sample these in our bootstrap calcu- 
lations this assumption will hold (so long as each of the sectors of 
sky is large enough to be representative). For this reason, in our sub- 
sequent analysis, we divide the SGP region of the 2dFGRS survey 
into eight sectors and the NGP into six. The selection of the regions 
has been made to ensure a statistically significant and roughly equal 
number of galaxies in each sector. These regions are then selected 
at random, with replacement, as in the standard bootstrap analysis. 
We make use of 20 bootstrap realisations in the analyses that fol- 
low, and use these to estimate the covariance matrix for each of our 
fits. The limitation of the bootstrap approach is that the samples are 
drawn from the observed volume of space, and may not represent 
the entire cosmic scatter. However, we find that error bars on pa- 
rameters derived from the mocks using the procedure of Hawkins 
et al. (2003) are in reasonable agreement with those derived from 
our bootstrap approach (see Table 2). 

3.2 The real-space correlation function 

Because the various redshift distortions to the correlation function 
only affect its measurement along the line-of-sight, it is possible to 
make an estimate of the real-space correlation function ^(r) by first 
projecting the two-dimensional correlation function, £(<r, n), onto 
the 7r = axis. This projected correlation function, H(cr), is given 

by, 

(• 1X3 

H(a)=2/ £(o,n)d*. (7) 
J o 

In practice the upper limit of the integration is taken to be a large 
finite separation for which the integral is found to converge. We 
find that limiting the integration to tt = 70 /i _1 Mpc suffices for 
the analysis presented here, providing us with stable projections 
out to a = 30 h' 1 Mpc. 

The projected correlation function can then be written as an 
integral over the real-space correlation function (Davis & Peebles 
1983), 

3(a) = 2 r rqr)dr _ (g) 

If we assume a power law form; £(r) = (V/ro) -7 , we can solve 
this equation for the unknown parameters, 




O0.1 1 10 



Projected separation ct [h Mpc] 

Figure 2. The projected correlation function, E(<x) /a, is shown for both 
relatively passive and active galaxies in the 2dFGRS. It can be seen that the 
correlation function of both sets of galaxies has an approximate power law 
form for a large range of separations, and this is illustrated with the best- 
fitting power law determined from this data (solid line). The dashed lines 
are extrapolations of these fits to larger and smaller scales. 



The projected correlation functions, H(a)/cr, are shown in 
Fig- HI together with error bars derived from the bootstrap reali- 
sations. It can be seen that for both sets of galaxies the power law 
assumption we have made is justified on small scales, and is quite 
consistent with the observations on large scales. The results from 
fitting a power law form for £ (r) are given in TableQ (Method P), 
together with those derived for the combined sample of all spec- 
trally typed galaxies. Upon comparing with Hawkins et al. (2003), 
we find that our estimate of the combined correlation function is 
essentially identical - from which we can conclude that restrict- 
ing our analysis to only those galaxies with spectral types has not 
biased our results in any noticeable way. 

In order to independently verify our assumption of a power 
law £(r), we have also calculated the real space correlation 
function, using the non-parametric method of Saunders, Rowan- 
Robinson & Lawrence (1992) (Fig. It can be seen that this 
method also estimates a power-law form for out to scales 
of ~ 20 h~ x Mpc and the best-fit parameters are shown in Ta- 
ble Q (Method /). We note however that a clear shoulder ap- 
pears to be present in both the correlation functions for separations 
r ~ 8 h~ x Mpc. It has recently been suggested that this may reflect 
the transition scale between a regime dominated by galaxy pairs in 
the same halo and a regime dominated by pairs in separate halos 
(e.g. Zehavi et al. 2003; Magliochetti & Porciani 2003). 

To increase the accuracy of our results we have split the galaxy 
sample into only two sub-samples. To justify this choice we have 
repeated our analysis on two further sub-samples of the active sam- 
ple. We found that the £(r) estimates were essentially identical (and 
consistent within the estimated uncertainties), and so the cluster- 
ing statistics are relatively insensitive to the exact amount of star- 
formation occurring in active star-forming galaxies. 
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Figure 3. The non-parametric estimates of the real-space correlation func- 
tions are shown for both our spectral types, using the method of Saunders, 
Rowan-Robinson & Lawrence (1992). It can be seen that our assumption of 
a power law form for £(r) is justified out to scales of up to 20 h^ 1 Mpc. 
The solid lines are the best-fitting power law fits shown in TablelTI whereas 
the dashed lines are extrapolations of these fits. 



Table 1. The derived parameters of the projected real-space correlation 
function. The fits have been determined using the range of separations 
0.2 < r < 20 h~ 1 Mpc. The Method column refers to the projected (P) 
or inverted (/) values of the parameters. Also shown are the values of o"g L 
derived from these correlation functions. Note that the lower limit of the 
fits (0.2 h^ 1 Mpc) was imposed to avoid biases from fibre collisions (see 
Hawkins et al. 2003). 



Galaxy type 


Method 


rn (h 1 Mpc) 


7 


NL 


All 


P 


4.69 ±0.22 


1.73 ±0.03 


0.83 ± 0.06 


Passive 


P 


6.10 ±0.34 


1.95 ±0.03 


1.12 ± 0.10 


Active 


P 


3.67 ±0.30 


1.60 ±0.04 


0.68 ± 0.10 


All 

Passive 
Active 


I 
1 
1 


5.01 ±0.23 
5.97 ±0.29 
4.12 ±0.32 


1.64 ±0.03 
1.93 ±0.03 
1.50 ±0.04 


0.88 ±0.05 
1.09 ± 0.08 
0.75 ± 0.09 



3.3 Relative bias 

The term bias is used to describe the fact that it is possible for 
the distribution of galaxies to not trace the underlying mass den- 
sity distribution precisely. The existence of such an effect would be 
a natural consequence if galaxy formation was enhanced, for ex- 
ample, in dense environments. The simplest model commonly as- 
sumed (somewhat ad-hoc) to quantify the degree of biasing present, 
is that of the linear bias parameter, b, 



= b\ 5 -P 



galaxies 



o 01 



0.1 



1 10 
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Figure 4. The relative bias between the most passive and actively star- 
forming galaxies is shown, in terms of (the square-root of) the ratio of the 
real-space correlation functions of these two samples. 



where p is a measure of the density, of either the mass or the galax- 
ies. A more specific model, based on the statistics of peaks (Kaiser 
1984; Bardeen et al. 1986), is that the degree of clustering we ob- 
serve in our galaxy sample, quantified in terms of the correlation 
function, £(r), is related to the mass correlation function in terms 
of, 



^(^)galaxies 



6 2 £(r) 



(11) 



Where here b is a constant that does not vary with scale, but more 
generally may depend on r. 

It is possible to estimate the magnitude of the biasing present 
in a sample of galaxies through the use of 'redshift-space distor- 
tions', and this issue will be addressed later in this paper. However, 
before proceeding it is already possible for us to determine the de- 
gree of relative biasing between our galaxy types at different scales, 
since, 

^passive (0 £passive(r) 



(10) 



- < ~TT ' (12) 

Active M 4activc(r) 

This relative bias between our two samples is shown in Fig. [4] 
where we have taken the ratio between the two estimates of the real- 
space correlation function, derived in the previous section. It can be 
seen that on small scales the clustering of the most passive galaxies 
in our sample is significantly larger than that of the more actively 
star-forming galaxies. The relative bias then appears to decrease 
significantly until on scales greater than about ~ 10 h^ 1 Mpc both 
samples display essentially the same degree of clustering (within 
the stated uncertainties). 

Another frequently used method of quantifying the degree of 
biasing present in a galaxy sample is through the parameter <jg L 
- the dimensionless standard deviation of (in this case) counts of 
galaxies in spheres of 8 h~ l Mpc radius. This quantity was deemed 
particularly useful because of the recognition (Peebles 1980) that 
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1, making the interpretation 



(13) 



for optically selected galaxies 0g' L 
of b particularly simple since, 

^ _ erg L (galaxies) 
o-g L (mass) 

Note that we write explicitly <t^ l to emphasise that this is a quan- 
tity defined on the nonlinear density field. It is an unfortunate stan- 
dard convention that, in the context of CDM models, as is used to 
denote an amplitude calculated according to linear theory. This is 
not the quantity considered here. 

With b defined in this way, the relative bias between our two 
spectral types is, 



frpassivc _ crg L (passive) 



(14) 



bactivc cr^ 1 - (active) 

The quantity o-g JL can be directly derived from our measured cor- 
relation functions in quite a straight-forward manner, since the ex- 
pected variance of the galaxy counts in a randomly placed sphere 
is. 



((N - Nf) R = N + 



(15) 



The first term is the shot-noise contribution and depends on how 
sparsely we have sampled the galaxy distribution. If we assume 
a power law form for the real-space correlation function we can 
estimate the fluctuation amplitude with the shot noise removed as, 



/ NL\2 



where, 



Mi) 



Mi) 



72 



(16) 



(17) 



[(3-7)(4- 7 )(6- 7 )27] ' 
(Peebles 1980). 

Our derived values of <7g L , for each of the galaxy samples 
considered in the previous Section, are given in Table □ It can 
be seen from this table that the relative bias of passive with re- 
spect to active galaxies (integrated over scales up to 8 h^ 1 Mpc) is 

bpassivo/kactivc = 1.09/0.75 = 1.45 ± 0.14. 

3.4 Modelling £(er, it) 

There is much further information to be derived from the observed 
£(<r, 7t) grids (Fig. 1) for each galaxy type. However, in order to 
do so we must first assume some model for the clustering of galax- 
ies with which to contrast the observations. Here we follow the 
analysis presented in Hawkins et al. (2003) with only minor modi- 
fications. Because the most significant limitation to this analysis is 
inevitably the assumptions that must be enforced upon our model, 
we summarise here the most important aspects of this model. 

In order to derive a model to fit the observed £(<r, n) grid, we 
need three main ingredients. The first is to assume some form for 
the real-space correlation function Because we are only go- 
ing to be concerned with relatively small-scale separations between 
galaxies ($; 20 h~ Mpc), we shall assume a power-law form of 
this function, 



£(r) = 



(18) 



In converting from real space to redshift space the next step is 
to account for the distortions in the correlation function which are 
caused by the linear coherent in-fall of galaxies into cluster over- 
densities (Kaiser 1987; Hamilton 1992), combined with non-linear 



velocity dispersion (e.g. Peacock et al. 2001). The linear-theory in- 
fall distortion can be written as (Hamilton 1992): 



C'(<7, 7T) = £o(*)Po(p) + Us)P2(p) + US)P^) 



(19) 



where Pt(fl) are Legendre polynomials, /i = cos(#) and 6 is the 
angle between r and n. The relations between £(r) and j3 for a 
simple power-law £(r) = (r-/ro)~ 7 are 

= (l + f + f) CM , (20) 



Us) = 



U(s) 



4/? 2 
7 



8(3 2 f 7(2 + 7) 



35 V(3- 7 )(5-7) 



CM 



(21) 



(22) 



We use these relations to create a model f ' (a, tt) which we then 
convolve with the distribution function of random pairwise mo- 
tions, f(v), to give the final model (Peebles 1980): 



£'(0-,7T- v/H )f(v)dv 



(23) 



and we choose to represent the random motions by an exponential 
form, 

where a is the pairwise peculiar velocity dispersion (often known as 
(J12). An exponential form for the random motions has been found 
to fit the observed data better than other functional forms (e.g. Rat- 
cliffe et al. 1998; Landy 2002). 

The factor /3 = n° m e /b arises from the growth rate in linear 
theory, 

f=**i„CS* t (25) 
a In a 

which is almost independent of the cosmological constant (Lahav 
et al. 1991), combined with the scale-independent biasing parame- 
ter b. A simple consequence of this model is that the redshift-space 
power spectrum will also appear to be amplified compared to its 
real-space counterpart 1 . It is worth explicitly re-stating that all of 
these derivations are based upon the assumptions of the linear the- 
ory of perturbations and linear bias, and assume the far-field ap- 
proximation (although this is not a significant issue for the 2dF- 
GRS, for which 

•^median — 0.1). 
It is also interesting to note that there is another cosmologi- 
cal effect that can result in the flattening of the observed £(cr, n) 
contours. It was first noted by Alcock & Paczynski (1979) that the 
presence of a significant cosmological constant, A, would result 
in geometric distortions of the inferred clustering if an incorrect 
geometry was assumed. However, Ballinger, Peacock & Heavens 
(1996) have shown that this is likely to be negligible for the low 
redshift data-set being considered here, for our assumed model of 

n m = 1 - n A = 0.3. 

To summarise, the parameters of our model are therefore the 

1 More precisely, the redshift-space distortion factor, /3, depends on the 
auto power spectra P mm (k) and P gg (k) for the mass and the galaxies, and 
on the mass-galaxies cross power spectrum P mg (fc) (Dekel & Lahav 1999; 
Pen 1998; Tegmark et al. 2001). The model presented here is only valid 
for a scale-independent bias factor b that obeys P gg (k) = bP rng (k) = 

b 2 Pmm(k). 
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real-space correlation function, £(r) = (r/ro) /3, and f(v), 
parameterised in terms of the velocity dispersion a. The best-fitting 
parameters of this model can now be determined from the model 
correlation function that best matches the observed £(er, n). 



4 RESULTS 

4.1 Validation of assumptions 

We have calculated the real-space correlation function, £(r), inde- 
pendently using the non-parametric method of Saunders, Rowan- 
Robinson & Lawrence (1992), in order to confirm the range over 
which it is a power-law (see Fig.[3J. For both samples of galax- 
ies, £(r) is adequately fit by a power-law to separations of r < 
20 ft -1 Mpc. This limit provides the upper bound to which we can 
compare the observed £(<r, it) with our assumed model. In addi- 
tion, because we have assumed the linear theory of perturbations in 
deriving our model we must impose a lower limit to the separations 
that we will use in our fit. To ensure that we have a sufficiently 
large fitting range we set this lower-limit to s — 8 h^ 1 Mpc. In 
fact it is quite plausible that the assumptions of linear theory are 
no longer valid at this separation. However, as will be shown, the 
ability of our model to recover the observed £ (a, 7r) at these scales 
is reassuring. 

On large scales it is known that the correlation function must 
deviate from a pure power law form, and there is some evidence 
to support this in Section lX2l on scales ~ 20 ft -1 Mpc. A num- 
ber of methods to account for this expected curvature in £(r) were 
investigated in Hawkins et al. (2003). However, the analysis pre- 
sented there suggests that so long as we restrict our fitting range to 
r < 20 h~ 1 Mpc the curvature has a negligible impact upon our 
parameter estimation. For this reason we neglect the possibility of 
curvature in the present analysis and restrict ourselves to using the 
simple power-law form for f (r). 

The other major assumption we have made is that the peculiar 
velocity distribution, f(v), has an exponential form. This can be 
tested using the method outlined by Landy, Szalay & Broadhurst 
(1998). This method makes a non-parametric estimate of the veloc- 
ity distribution, using the Fourier decompositions of the observed 
f (a, 7r) grid along the k a = and = axes. Unfortunately this 
method ignores the effects of coherent in-fall, which can substan- 
tially change the resulting estimate of a (see Hawkins et al. 2003). 
However, it is found that the recovered f(v) is well fit by an ex- 
ponential form - for both types of galaxies. Hawkins et al. (2003) 
have shown that although incorporating the effects of coherent in- 
fall changes the estimated velocity dispersion, a, substantially, the 
method gives a robust estimate of the form for /(«). 

4.2 Parameter fits 

All four of our parameters (/3, a, ro and 7) are allowed to vary 
over a large range of possibilities and a downhill simplex multi- 
dimensional minimisation routine is adopted to find their best- 
fitting values (see e.g. Press et al. 1992). Our calculated £(cr, n) 
contours are shown in Fig. [5] together with those of the best-fitting 
model correlation functions derived in this manner. The peak pa- 
rameters of this best-fitting model are detailed in Table. |2| together 
with their estimated uncertainties. We find that there is quite a sig- 
nificant degeneracy between (3 and a (see Fig.|6j. This is also exac- 
erbated by the relatively noisy nature of £(r) at these scales, which 
makes ro and 7 difficult to constrain accurately. 
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Figure 6. The estimates of /3 and a for the bootstrap samples are shown 
for both the passive (squares) and active (Mangles) galaxy samples. It can 
be seen that a significant degeneracy exists between (3 and a, for both 
samples. The fits have been made using only the quasi-lineat regime of 
8 — 20 /i _1 Mpc. The crosses show the fits to the full samples together with 
the 1(7 uncertainties shown in Table 2. 



One immediate conclusion is that the velocity dispersions of 
the two galaxy populations are very distinct, even taking into ac- 
count the substantial statistical uncertainties. This is an interesting 
result which has significant implications for the proportion of each 
of these galaxy types we expect to occupy large, virialised clusters 
of galaxies. 

Another conclusion that we can easily make is to quantify the 
relative bias between our two spectral types, as described in Sec- 
tion !3.3l However, as demonstrated in that Section, the relative bias 
between our galaxy types is in fact essentially unity over the range 
for which our model assumptions are valid (8 — 20 h^ 1 Mpc), a 
result confirmed in this analysis. However, a much more important 
quantity that can be inferred from these redshift-space distortions, 
that could not be determined previously, is the absolute value of the 
biasing between the galaxy and mass distributions, b, as described 
in Section l3~3l We return to this point in the next Section of this 
paper. 



5 COMPARISON WITH PREVIOUS RESULTS 

Because previous estimates of j3 have used slightly different galaxy 
samples, it is first necessary to correct for various effects before 
making a proper comparison (see Lahav et al. 2002). There are 
two main issues which affect the different estimates of the bias- 
ing; the effective redshift of the survey sample used, z s , and the 
effective luminosity of the galaxies in that sample, L s . These quan- 
tities vary between the samples used, depending on the weighting 
scheme adopted and the limiting redshift of the survey. 

Because we have only used 2dFGRS galaxies for which a 
spectral type is available our sample is limited to z max = 0.15. 
To determine the effective redshift of our sample it is necessary for 
us to determine the weighted average of the galaxies used in each 
of our calculations. Doing so reveals that for all three of our sam- 
ples z s = 0.11. In a similar way we can calculate the weighted 
mean luminosity of each of our samples, which are found to be as 
follows. Combined: L s = 1.06L,; Passive: L s — 1.26L,; Active: 



8 D.S. Madgwick et al. 



Table 2. The best fitting model parameters derived from the observed 7r) grid are shown. All fits have been made over the quasi-linear redshift-space 
separation range 8 < s < 20 h" 1 Mpc. The quoted uncertainties correspond to the \a scatter derived from the bootstrap estimates (see Fig.[S). For 
comparison, uncertainties in /3 have also been estimated from sparsely sampled 2dFGRS mock galaxy catalogues (Cole et al., 1998) limited to z < 0.15, and 
correspond to A/3 = 0.12 for the full sample and A/3 = 0.16 for a l-in-2 random sampling. This demonstrates that the bootstrap approach has given a fair 
assessment of the cosmic scatter in these estimates. 



Parameter 


All galaxies 


Passive galaxies 


Active galaxies 


P 


0.46 ± 0.10 


0.48 ±0.14 


0.49 ±0.13 


a 


537 ±87 km s" 1 


612 ±92 km s- 1 


416 ±76 km s- 1 


ro 


5.47 ±0.32 ft- 1 Mpc 


7.21 ± 0.34 h- 1 Mpc 


4.24 ± 0.41 h- 1 Mpc 


7 


1.75 ± 0.08 


1.91 ±0.10 


1.60 ±0.11 




Table 3. Comparison between biasing results derived using the 2dF Galaxy Redshift Survey by various authors. The results of Peacock et al. (2001) are derived 
from the redshift-space distortions in the two-point correlation function (Note that this galaxy sample was much deeper as galaxies without spectral types were 
used). Lahav et al. (2002) made their estimate of the bias through comparing the amplitude of fluctuations in both the 2dFGRS and the CMB. Verde et al. 
(2002) calculated the bi-spectrum of the 2dFGRS, which constrained the linear bias parameter, b, which we have converted to /3 by assuming our concordance 
cosmological model of a flat Universe with Qm = 0.3. Note that the results of Lahav et al. (2002) and Verde et al. (2002) are valid over scales expressed in 
terms of wavenumber k rather than real-space distance. We have converted between the two by simply taking r ~ 1/fe. 



Galaxy type 


Author 


Scales (h 


^Mpc) 


(z s , Ls/L*) 




13(0, L,) 


All 




8 - 


20 


(0.11,1.06) 


0.46 ±0.10 


0.44 ± 0.09 


Passive 




8 - 


20 


(0.11,1.26) 


0.48 ±0.14 


0.47 ±0.14 


Active 




8- 


20 


(0.11,0.95) 


0.49 ±0.13 


0.48 ±0.13 


All 


Hawkins et al. 


8 - 


20 


(0.15,1.4) 


0.49 ± 0.09 


0.47 ± 0.08 


All 


Peacock et al. 


8- 


25 


(0.17,1.9) 


0.43 ± 0.07 


0.45 ±0.07 


All 


Lahav et al. 


7- 


50 


(0.17,1.9) 


0.48 ± 0.06 


0.50 ±0.06 


All 


Verde et al. 


2 - 


10 


(0.17,1.9) 


0.56 ± 0.06 


0.59 ± 0.06 
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L s = 0.95L,; where we have taken M* — 51og 10 (/i) = —19.66 
(Norberg et al. 2002b). 

Assuming linear dynamics and linear biasing, the redshift- 
distortion parameter, f3, for a given sample redshift and luminosity 
can be written as, 



(3(L,z) 



b(L,z) 



(26) 



The evolution of the matter density parameter, f2 m (z) is straight- 
forward to determine, assuming a given cosmological model, 

Q m (z) = On(l + zf(H/H y 2 , (27) 
where fi m is the matter density at the present epoch and, 
/ H \ 2 

= + zf + (i - n m - n A )(i + zf + n A . (28) 

The determination of the variation in the biasing parameter, 
b, with redshift is much less straight-forward. As shown in Sec- 
tion !3.3l b can be defined as, 

-~ NL ' s (z) 



b(z) = 



r^(z) ' 



(29) 



where here we have added a redshift dependence, b(z), and labelled 
the two crg' L 's by the superscripts g and m to denote galaxies and 
mass respectively. As described by Lahav et al. (2002), there is now 
much evidence to suggest that whilst the matter fluctuations con- 
tinue to grow at low redshifts, the fluctuations in the galaxy distri- 
bution are relatively constant between < z < 0.5 (see e.g. Shep- 
herd et al. 2001). If we assume that the matter fluctuations grow 
according to the linear theory of perturbations then, Og' m (z) = 
(7§ Ll m (0).D(z), where D(z) is the growing mode of fluctuations 
(Peebles 1980). Whereas, cr£ L ' s (£,z) w al 



j^ L ' g (L,0). Therefore, 



b(L,z) 



b(L,0) 
D(z) 



(30) 



The final step then is to correct the biasing parameter, b, for the 
luminosity of our sample. Norberg et al. (2001) found from the 
analysis of the galaxy correlation functions on scales < 10 h^ 1 
Mpc that, 



0.85 + 0.15 



(3D 



b(L y 0) 
6(X*,0) 

Assuming that this relation also holds in our quasi-linear regime of 
8 — 20 h^ 1 Mpc, then allows us to determine f3 at redshift 2 = 
and luminosity L — L*. 2 

Table [3] shows the results for j3 derived in the analysis pre- 
sented here, both before and after converting to redshift z — and 
luminosity L — L t . Also shown are other results derived from the 
2dFGRS by previous authors. It can be seen that there is a remark- 
ably good agreement between all the results presented. We note that 
these results have been derived by applying linear corrections to a 
selection of quasi-linear regimes, which may introduce systematic 
errors into our results. This is a particular concern for the results 
of Verde et al. (2002), which correspond to the smallest separation 
ranges used. 

2 Note that in converting the linear bias parameter, b(L, 0), for each of 
our spectral types to b(L t , 0), we have explicitly assumed that each type 
of galaxy displays the same variations in clustering with luminosity. This 
result has been verified by Norberg et al. (2002a), who calculated the clus- 
tering amplitudes for different galaxy samples divided in spectral type and 
luminosity. 



6 DISCUSSION 

We have derived a variety of different parameterisations for the 
2dFGRS correlation function, £(<r, n), for different spectral types. 
The two types we have used can roughly be interpreted as divid- 
ing our galaxy sample on the basis of their relative amount of cur- 
rent star-formation activity, and hence provide useful insight into 
how galaxy formation may relate to the large-scale structure of the 
galaxy distribution. The actual cut we have imposed is most natu- 
rally interpreted in terms of the Scalo birthrate parameter, &s ca i . 
This is defined to be the ratio of the current star-formation rate and 
the past averaged star-formation rate. Adopting this convention, our 
cut of 77 = —1.4 corresponds to dividing our sample into galaxies 
with fcscaio = 0.1, i.e. between galaxies whose present star forma- 
tion rate is greater or less than 10% of their past averaged rate. 



6.1 Relative bias on small scales 

On scales smaller than ~ 8 h^ 1 Mpc the clustering of pas- 
sive galaxies is much stronger than that of the more actively star- 
forming galaxies. This was demonstrated quantitatively by the real- 
space correlation functions derived in Section l3~2l for which the 
passive galaxy sample were fit by a power-law with larger scale 
length, ro and steeper 7. In addition it was shown that the values 
of erg L derived for each of these samples were quite distinct, being 
o£ L = 1.09±0.08 for the passive galaxies and o£ L = 0.75±0.09 
for the actively star-forming galaxies, implying an (integrated) rel- 
ative bias between our two types of, 



= 1.45 ±0.14 , 



(32) 



at the effective redshift and luminosity of our galaxy samples (see 
Table 3). Note that this ratio quantifies the integrated relative bias 
between scales of — 8 h^ 1 Mpc. 

Our correlation functions per type confirm that the slope of 
passive (early type) galaxies is steeper than that of active (late tape) 
galaxies (cf. Zehavi et al. 2001). On the other hand, the slope of the 
correlation functions derived for different luminosity ranges show 
no significant variation (Norberg et al. 2001, 2002a; Zehavi et al. 
2001). These results call for theoretical explanations and they set 
important constraints on models for galaxy formation. 



6.2 Velocity dispersions 

The velocity distributions of our two samples were found to be dis- 
tinct. The passive galaxy sample displayed a consistently larger ve- 
locity dispersion, a, than the actively star-forming sample on all 
scales, and in particular on separations of 8 — 20 h~ x Mpc were 
found to be 612 ± 92 and 416 ± 76 km s^ 1 respectively. This result 
is consistent with the observations of Dressier (1980), that a signif- 
icant morphology-density relation exists - since a larger velocity 
dispersion would tend to suggest a higher proportion of galaxies 
occupying virialised (high-density) clusters. 



6.3 Relative bias on large scales 

The determination of the redshift-distortion parameter, j3, was 
found to be much less straight-forward. The evidence from our 
analysis is that j3 has only a relatively small dependence on the 
spectral type of the galaxy sample under investigation. We found 
that on scales of 8 — 20 h^ 1 Mpc, the two redshift distortion param- 
eters were; /3 = 0.48 ± 0.14 and f3 = 0.49 ± 0.13 for the passive 
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and actively star-forming galaxy samples respectively, yielding a 
relative bias of only, 



^active 

The overall redshift-distortion parameter, /3, independent of spec- 
tral type is found here to be 0.46 ± 0.10 (on scales of 8 — 
20 h~ x Mpc), at our sample's mean redshift and luminosity. By 
making various assumptions (Section|5J this result can be converted 
to redshift z = and L„ luminosity, giving /3(0) = 0.44 ± 0.09. 
This result is is almost identical to the /3(0) = 0.47 ± 0.08 derived 
from the results of Hawkins et al. (2003), using the entire 2dFGRS 
data-set, over the same separation range. 

In the analyses presented in this paper two fundamental limits 
were found to greatly inhibit our ability to accurately characterise 
the relative and absolute biases on different scales. The first of these 
was that on small scales - where the clustering of our two popula- 
tions is most distinct - the assumptions of our model of the galaxy 
clustering were no longer accurate, and so we could not accurately 
determine ft or a on these scales. Our second limitation was found 
on large scales (s ~ 20 /i -1 Mpc), where the galaxy correlation 
functions became noisy and were no longer well parameterised by 
a power-law form. 

The latter of these issues can be addressed to some degree 
simply by a change of formalism to incorporate the power spec- 
trum estimations of each galaxy type or colour (Peacock 2003). 
Because this characterisation of the clustering is more sensitive to 
larger scales of separations it would allow us to more rigorously 
test whether the large-scale (s > 20 h" 1 Mpc) clustering of these 
populations are in fact distinct and also allow us to incorporate the 
possibility of scale-dependent bias. The derived correlation func- 
tions per type could also be used within the framework of halo oc- 
cupation number to derive e.g. the mean number of galaxies of a 
given type per halo (e.g. Zehavi et al. 2003; Magliochetti & Por- 
ciani 2003). 
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